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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

INVEST 3GATION OF A 10-STAGE SUBSONIC AXIAL-FLOW RESEARCH COMPRESSOR 

I - AERODYNAMIC DESIGN 
By Irving A. Johnsen 


SUMMARY 

A 10-stage axial -flow compressor was designed for use as a research, 
unit in which the problems associated with the compounding of high per- 
formance stages could be studied. The design was based on a symmetrical 
velocity diagram with constant total enthalpy at all radii. Selected 
design limitations included a Mach number limit of 0.7, a loading limit 
(Cptf where Cp is the lift coefficient and a is the blade-element 

solidity) which increased progressively from 0.8 to 1.0 through the com- 
pressor, a solidity limit of 1.1, an inlet hub-tip ratio of 0.55, and a 
tip diameter of 20 inches. The application of these design limits with 
an assumed stage polytropic efficiency of 0.90 gave the following design 
values: a total-pressure ratio of 6.45 and a weight flow of 57.5 pounds 

per second at the design tip speed of 869 feet per second. 

The over-all performance characteristics of this compressor, as 
determined experimentally, are included. 


INTRODUCTION 

The desirability of compactness in axial -flow compressors for gas- 
turbine power plants requires that multistage units produce high pres- 
sure ratios per stage and high mass flows per unit frontal area, while 
at the same time maintaining high efficiencies. Two-dimensional high- 
speed cascade results (reference l) indicated that stage pressure ratios 
could be realized which were considerably in excess of those currently 
being obtained in multistage machines. In order to study the possibility 
of utilizing these high loadings and high relative Mach numbers, a rotor- 
blade row was designed and Investigated (reference 2). The results of 
this investigation verified the fact that Increased pressure ratios were 
attainable in a single stage without significant penalties in efficiency. 

Even with the realization of increased pressure ratio in a single 
stage, however, it was questionable whether the results were applicable 
to multistage units. Additional factors which are introduced by multi- 
staging, such as boundary- layer growth, oscillatory radial motion of the 
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gas through many stages, blade-row matching, off -design operating 
problems, and the like, may prevent high- stage performance from being 
realized. Therefore, a multistage axial -flow compressor was designed 
and constructed in which the attempt was made to achieve higher stage 
pressure ratios than those in current use, while at the same time main- 
taining high mass flow and efficiency. This design was intended to 
provide a research unit in which the problems arising from the com- 
pounding of high performance stages could be studied. Inasmuch as it 
was not feasible to attempt to attain an absolute max imum of perfor- 
mance in each stage of this initial design, design limitations were 
chosen such that there was reasonable assurance that the desired high- 
performance multistage research compressor would be obtained. 

The compressor design, including the deter minat ion of the velocity 
diagram and the selection of the compressor blading, was conducted in 
1947 at the HACA Lewis laboratory based on the best Information avail- 
able at that time. This report summarizes the design procedure and 
the design characteristics of this 10-stage compressor unit. The over- 
all performance map of the compressor for the operating range from 
30 percent to 110 percent of equivalent design speed is included. 


SYMBOLS 

The following symbols are vised in this report: 

A area (sq. ft) 

a local velocity of sound (ft/sec) 

b intercept in equation Ap = mx + b 

C L lift coefficient 

c L,i design camber (theoretical Cl for Isolated airfoil) 

Cx,cr loading limit 

g acceleration of gravity (ft/sec^) 

h dimensionless ratio of axial component of air velocity to rotor- 

tip speed 

h average dimensionless axial velocity across passage 

hp horsepower 

K constant 
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M average axial Mach number across passage 

m slope in equation Af3 = m + b 

n polytropic exponent 

P total or stagnation pressure (lb/sq ft) 

p static or stream pressure (lb/sq. ft) 

Q volume rate of flow (cu ft/ sec) 
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V 2 

q. velocity head, (ft) 

2g 

R gas constant 

r radius (ft) 

T total or stagnation temperature (9R) 

t static or stream temperature (°R) 

U velocity of rotor blade element at radius r (ft/sec) 

V absolute air velocity (ft/ sec) 

W air velocity relative to rotor (ft/sec) 

Wf air weight flow (lb/ sec) 


X ratio of circumferential component of absolute velocity of 

entering air to rotor -tip speed 

Y ratio of change in circumferential component of velocity to 

rotor-tip speed 

Z ratio of blade-element velocity to rotor-tip speed (compressor 

radius ratio) 

a angle of attack (deg) 

{3 absolute inlet-air angle, measured with respect to axis (deg) 

(fig. 3) 

inlet -air angle relative to rotor, measured with respect to axis 
(deg) (fig. 3) 


P' 
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r 

Ap 

n 

\ 

p 

p t 

a 

Subscripts: 

e 

h 


ratio of specific heats, Cp/cy 
turning angle (deg) 
adiabatic efficiency 
polytropic efficiency 

static or stream density (lb/cu ft) 
total or stagnation density (lb/cu ft) 

blade-element solidity, chord/ spacing 

effective value in velocity diagram corrected to con- 
stant social velocity 

hub 


m mean 

t tip 

0 circumferential direction 

0 station ahead of inlet -guide vanes 

1.3.5 ... 19 station ahead of rotors of 1 st , 2 n<3 - . 

2.4.6 . . . 20 station ahead of stators of l s ^, 2 ni 

21 station ahead of exit vanes 

22 station after exit vanes 


lO^ 1 stages 
. 10 th stages 


GENERAL CHARACTERISTICS AND LIMITATIONS 

Compressor size . - The attainment of a compressor of minimum size 
for a given pressure ratio and weight flow requires the consideration of 
velocity distribution as well as limitations on the turning of the air 
and the blade relative Mach number. The weight flow desired for this 
compressor was of the order of 50 pounds per second. On the basis of 
preliminary design estimates and with the use of the established aero- 
dynamic limitations, the inlet diameter was fixed at 20 inches. This 
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preliminary est ima te also established the fact that the desired pres- 
sure ratio of 6 could be obtained in 10 stages. 

Compressor configuration . - In a multistage axial -flow compressor, 
configurations are possible in which the tip diameter, the hub diameter, 
or both vary from stage to stage through the compressor. As pointed 
out in reference 3, a progressive increase in tip diameter is desirable 
for obta ining maximum pressure ratio per stage. However, this advan- 
tage is .balanced by the fact that (for a given inlet stage), an increase 
in tip diameter results in very short blades, high ratios of hub-to-tip 
diameter, and small aspect ratios for a given blade chord. Conse- 
quently, tip clearance and annulus losses would be larger than for a 
comparable constant tip-diameter design. Since construction would also 
be simplified, a constant tip diameter was selected. 

In order to maintain the limiting relative Mach number on each 
blade row (and thereby obtain a max imum pressure ratio), the hub 
diameter was progressively increased through the compressor. This 
served to increase the blade -element velocity at the hub, as well as 
the axial velocity of the stream, and therefore to maintain the limiting 
relative Mach number despite the increasing speed of sound through the 
unit. 


Velocity diagram . - The type of velocity diagram used in axial- 
flow compressor design is vitally important in that it establishes the 
blade element and stage performance. On the basis of a blade -element 
analysis, the optimum velocity diagram was found to be symmetrical 
(reference 3). 

In applying the blade-element theory to stage design, it was con- 
sidered desirable to maintain a constant total enthalpy from hub to 
tip in order to avoid an energy unbalance that might result in exces- 
sive mixing losses. It is shown in reference 3 that the use of a sym- 
metrical velocity diagram with constant total ent halp y gives M ghPT* 
pressure ratio and specific mass flow than the more conventional free- 
vortex design. This high performance results from the fact that the 
inlet Mach number is maintained near the imposed limit 'for all blade 
elements in both rotor and stator. 

In view of the generally desirable pressure ratio and mass-flow 
characteristics, the symmetrical velocity diagram with constant total 
enthalpy at all radii was selected for this design. 


t-wrrmrff .-ss***. 
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Stage H notations . - The stage performance of axial -flow compres- 
sors is dependent on the limitations imposed on the relative Mach 
n umb er of the air and the turning done by the blade rows. The Mach 
number limitation was applied by using a value of 0.7 for flow relative 
to the compressor rotor (the maximum value exists at the hub in this 
design). In the calculation used, this results in a Mach number at the 
following stator hub that is slightly higher than 0.7. This Mach num- 
ber level was chosen on the basis of performance results obtained with 
the HACA 8-stage compressor (reference 4.). Subsequent data (refer- 
ence 2) have shown that higher Mach numbers could have been used without 
serious losses in efficiency; however, the selected Mach number level 
provides a reasonably high pressure rise and insures high efficiency. 

The solidity at the hub for all blade rows was selected as 1.1. 

This value is of the order of that successfully used in the design 
reported in reference 4. 

There was little information available on the amount of turning 
that could be imparted to the air under various operating conditions. 

Therefore, a limitation on Cjcs (as was used in reference 4) was 

adopted because of convenience in calculations. In this particular ,, 

design, this limitation occurs at the hub and is essentially the 
same In rotor and stator. The values of C^cr were selected to increase 
progressively through the compressor (Cxp of 0.8 for the first three w 

stages, 0.9 for the next three stages, and 1.0 for the last four stages). 

This variation in the loading limitation through the compressor was 
selected to improve the off -design operation of the compressor, particu- 
larly at speeds lower than design and at low flow. Under these opera- 
ting conditions, the angle of attack in the first stages is increased. 

Since lightly loaded blades can be expected to have a wide operating 
range, the use of a low design value of C^a results in blades that 

will accept this increased angle of attack without flow breakdown. In 
addition, the change in density through the compressor at low speed is 
proportionally less than the design decrease in area; this successively 
decreases the angle of attack and therefore unloads the latter stages. 

These two effects result in more nearly uniform loading through this 
compressor at low speeds. A more recent and complete analysis of the 
off -design operating problem In multistage axial-flow compressors has 
Indicated that the optimum loading variation is probably one in which 
the loading increases in the middle stages and decreases again in the 
latter stages. However, this selected increase in C^c through the 

compressor was expected to provide some improvement in the characteris- » 

tic under low-speed and low-flow operation. 

Mass flow and pressure ratio for first stage . - For gas-turbine * 

power-plant application, a sms, 41 : • hup~diaineter- and high axial air velo- 
city are desirable to obtain as high a mass flow as possible consistent 
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with high pressure ratio. For convenience in establishing the design 
conditions for the inlet stage, the optimum combination of weight flow 
and pressure ratio was considered on the basis of useful power input to 
the air. In order to establish these conditions, the following power 
parameter (appendix A) was used: 


K(hp) = 1 

^t a l 



1.0 


For given inlet conditions and compressor diameter, and for given limi- 
tations of Mach number, lift coefficient, and solidity, this power 
parameter is a function of the hub-tip ratio and the dimensionless 

axial velocity of the air at the hub hjj. 

In figure l(a), the power parameter (determined on the basis of the 
limitations chosen to apply to the first stage) is shown as a function 
of Zh for various values of 1%. It can be seen that the thr-x-t - mum 
power is obtained at a hub -tip ratio of the order of 0.55 and that the 
curves are essentially flat from 0.50 to 0.55. The hub-tip ratio for 
the inlet stage of the compressor was selected as 0.55. 

As shown in figure l(b), maximum power for the hub-tip ratio of 
0.55 is obtained at a dimensionless axial velocity hfr of approximately 
0.70. A further consideration that must enter into the choice of h^ 

at the inlet, however, is the axial -velocity distribution over the rest 
of the stage. A disadvantage of this symmetrical diagram, constant 
total-enthalpy design, is the fact that the axial velocity at the rotor 
tip, particularly at the exit, may become dangerously low. The axial 
velocity of the air discharged at the rotor tip h2,t was therefore 
investigated by use of the simplified radial equilibrium approximation, 
which considers only the balance of centrifugal and pressure forces and 
neglects the effect of curvature of the streamlines. In order to keep 
the axial velocity at the rotor tip at a reasonable level, a value of 
hfc. of 0.75 was selected for this design. This provides a higher mass 
flow and a power Input only slightly less than the maximum attainable. 

Efficiency . - As pointed out in reference 3, the exact analysis of 
efficiency of an axial -flow compressor is extremely complex and must be 
based upon a detailed knowledge of the flow processes involved. How- 
ever, an indication of the efficiency of the main portion of the flow 
can be obtained by a consideration of blade -element efficiency. An 
examination of the velocity diagrams on the basis of blade -profile 
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losses (fig. 17 of reference 4) showed that the blades of this selected 
entrance stage operate in the region of - peak efficiency with the lowest 
blade -element efficiency existing at the hub. 

A preliminary estimate was made to evaluate the efficiency of 
succeeding stages on a blade-element basis. This study showed that the 
shift in blade -element efficiency in succeeding stages is in the direc- 
tion of increased efficiency. It was recognized that viscous effects 
probably become more severe in later stages and that these real flow 
effects modify the blade -element trend. However, on the basis of a 
blade-element consideration of pressure ratio, weight flow, and effi- 
ciency, the choice of velocity distribution appeared reasonable. 

Axial- velocity distribution . - In this design, which is based on a 
symmetrical velocity diagram with constant total enthalpy and simplified 
radial equilibrium, axial velocities are greatest at the hub. This 
characteristic permits the use of large values of TCJt over the entire 
blade height without exceeding the Cjp limi tation, and allows a Mach 
number near the limiting value to occur across the entire blade height, 
both of which contribute to high performance. The disadvantage of this 
design, however, is that large axial -velocity changes may take place 
across the blade rows. Across the rotor tip, for example, a decelera- 
tion occurs. If this deceleration is sufficiently rapid, flow separa- 
tion on the blade surface may result. 

In order to establish design velocity diagrams, axial -velocity 
changes must be evaluated . It is shown theoretically in reference 5 
that the effect of radial motion of the gas may be important in the 
determination of this axial -velocity distribution, particularly when 
the blade aspect ratio is large. However, the exact determination of 
axial -velocity distribution (considering the complete radial equilibrium 
equation) requires a long process of step-by-step calculation, which is 
somewhat impractical in the case of a multistage compressor. Therefore, 
for this design, the effect of radial motion was neglected and the axial- 
velocity distribution (at the stations imme diately upstream and down- 
stream of blade rows) was determined solely on the basis of simple equi- 
librium of pressure and centrifugal force. 

Subsequent investigation (reference 6) indicates that the analysis 
based on simple radial equilibrium gives a reasonable approximation to 
actual conditions across the rotor. Furthermore, analysis has shown 
that across a complete stage the differences between the complete radial 
flow solution of reference 5 and the simple radial equilibrium approxi- 
mation are small. It therefore appears that the use of simple radial 
equilibrium relations provided a sufficiently accurate determination of 
the axial -velocity distribution. The simple radial equilibrium rela- 
tions ere developed in appendix B. Charts for determining average 
dimensionless axial velocities on this basis are presented in figure 2. 
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Average axial component of velocity . - In high -performance axial- 
flov compressors, the problem of reducing the velocity to the value that 
is required for combustion is a difficult one. A possible solution is 
to reduce the average axial velocity in passing through the compressor; 
this would reduce the requirements of the diffuser. It would also 
increase the blade height (particularly in the latter stages), which 
would be desirable because of reducing tip clearance and annulus losses. 
However, this method has three disadvantages: 

(1) The relative Mach number is reduced, thus reducing the stage 
pressure ratio and requiring more stages for a given over-all pressure 
ratio . 

(2) The increase in blade height results in a more extreme axial- 
velocity profile across the passage of the latter stages, which might 
result in increased mixing losses in the process of being converted to 
a uniform flow in the diffuser. 

(3) The value of Ap/q across a given blade row will be increased, 
which would have an adverse effect with regard to boundary-layer build- 
up. 


In view of these disadvantages, the design was carried out to main- 
tain the relative Mach number in all stages with one limitation. The 
average axial Mach number was not permitted to exceed 0.6. 


DESIGN CALCULATION METHOD 

The basic theory used in the velocity diagram determination Is 
essentially a method proposed by J. Austin King, with velocities 
expressed nondimensionally as ratios of the rotor tip speed. A typical 
velocity diagram for the case where axial velocity is constant is shown 
in figure 3(a). In the design calculation, however, the axial velocity 
varies through the stage because of (l) variation in hub diam eter, 

(2) compressibility, and (3) the change in axial velocity required to 
satisfy the condition of simplified radial equilibrium between each 
blade row. A rotor velocity diagram for the case of a varying axi 
velocity is shown in figure 3(b). This diagram is corrected to con- 
stant axial velocity in order to permit use of cascade data; 1%, W e , 

Pe ' > Ag e are "effective" values. This procedure for correcting to 

an equivalent constant axial -velocity diagram is in accordance with the 
recommendation of reference 7. Reference 6 and unpublished data 
obtained at the Langley laboratory indicate that this correction system 
is adequate for most applications, although it may over correct when 
axial -velocity changes are large. 


i 
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Velocity diagram calculations . - The calculation procedure for 
determining velocity diagrams through the compressor waB as follows: 

Step 1: 

= 0.55 and 3% = 0.75 at the inlet of the first stage 

Step 2: 


C^cr = 0.8 for the first stage 


Step 3: 

h was selected from the design chart for the given values 

of Zfa, 1%, and Cpcr 


Step 4: 


Yh 




+ 



as derived in reference 4 


Step 5: 


Step 6: 


Step 7: 


Step 8: 


Y t = 


= 


*31 ' Y h 
2 


*t - 


1 


- Yt 
2 


^t _ Mach number limit _ 

a h " A/h h 2 + (Z* - X h ) 2 ' ^ + (Z* - X h ) 2 


^t 3^h 

&t ~ a h ®t 


Step 9: 


i 
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Substituting the relation derived in appendix B for a^/a-t gives 


u t 


1 

*t " 

M 

1 + rzi 

(Ut/a h ) 2 2 

1 — Zh 2 , p. 1 

[_ 4 , (1 + V) - *t lc >Se ZjJ 


Step 10 i 

The polytrqpic efficiency of this compressor was estimated as 0.90} for 
convenience in calculations, the polytropic exponent n was used where 



Step 3JL: 

Bernoulli's equation written in terms of static-pressure ratio is 

n 

n - 1 


An analysis was made on the effect of the velocity components X and h 
on stage pressure ratio, and it was established that 

x l 2 + h l Z * X 3 2 + *3 2 

This was also found true for succee ding stages; therefore, these terms 
could be deleted with little error. The pressure ratio was therefore 
determined as 


P3 

Pi 


1 + (r - 1) 



+ £l*2 


+ YZ . 


Step 12: 




P3 _ 
Pi " 


1.0 + (r 




n 


n - 1 




1 


n - 1 


YZi 


1.0 + (r - 1) 
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Step 13: 

^3 

t l 

Step 14: 

IJt = Ut/ai 
a 3 ^sAl 

The static-pressure ratio as determined in step 11 (as well as the 
density and temperature ratios) varied only slightly from hub to tip 
because of the variation in a]_, which was small. ■ Therefore, for design 

purposes, including the determination of the power parameter (appendix A), 
the pressure ratio was taken as that calculated at the tip station. 

Step 15 j 

At this point, a trial-and-error solution was required to select 
Z 3 h in order that both the Mach number limit at the next stage and 

the continuity equation were satisfied. 


/tr+N 

l.o + (r - i) Ur 


YZ 


_ __(l - z, h 2 ) 

h.3 = hi 1 — L 


Pl 


(l - Zt,h 2 ) \ P 3/ 


Step 16: 

Then, from the design chart (fig. 2(a)) for the assumed value of 
the value of 113 ^ was determined. 

Step 17: 

Steps (l) to (7) were repeated to satisfy the equation under 
step 8 at station 3. 

^3,h 2 + ( Z 3,h - X 3,h ) 2 (aj—) = °’ 7 

The calculation procedure then defined the tentative hub diameter, 
stage by stage, through the compressor. The calculations were carried 
through until the number of stages was sufficient to produce an over- 
all total-pressure ratio of approximately 6 . 
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With, the tentative hub diameters determined, axial chord lengths 
and gaps were selected and a curve of rotor hub diameters was established. 
This curve was then faired to eliminate profile discontinuities and to 
establish the actual diameters, for rotor and stator, through the com- 
pressor. The preceding calculations were then repeated for the new 
values of hub d iam eter in order to establish the final design values. 

This fairing of the hub diameters resulted in a small deviation 
from the 0.7 Mach number limit desired for the hub of the rotor. Act ual 
inlet Mach number values at other stations were established from the 
relation at the rotor inl et 

W _tf_ ^t 

a ' Ut a 

and, at the stator inlet 

V = JT Ufc 
a Ut a 

where, for the general case, the equation under step 9 can be written 
for rotor inlet 



and, for stator inlet 



Also, for the rotor inlet 


_W_ 

Pt 


2 


= (X + Y) 2 + h 2 
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and for the stator inlet 


( t ^) 2 " ( x + Y > 2 + * 2 

Elis distribution of inlet Mach, number across the annulus Indicated how 
well any particular blade row was being utilized. 

In selecting the hub diameter for a desired flow variation through 
the compressor, it is desirable to account for the presence of boundary 
layer. This, however, requires a precise knowledge of the magnitude of 
the boundary layer and its variation through the compressor. Because 
of a lack of pertinent data on the subject at the time of the design, 
no attempt was made to provide a boundary-layer allowance. It was 
therefore anticipated that the design weight flow might not be obtained. 
Furthermore, desired angles of attack might not be maintained, particu- 
larly in the latter stages where boundary layer (as a percentage of pas- 
sage height) is a maximum. However, an observed tendency to obtain 
greater rotor turning (and therefore pressure ratio) than predicted 
from low-speed cascade data would tend to provide a balancing effect 
with regard to angles of attack. 

Construction of velocity diagrams . - At any given radius, the rotor - 
inlet velocity vectors could be established from the data of the calcu- 
lation procedure and the radial equilibrium relations at the rotor inlet 
station. Eotor exit vectors for this same radius were then established 
from: (l) values of h2 (for the new value of Zh at the rotor exit) 
determined on the basis of continuity and simple radial equilibrium at 
this station and (2) values of X and Y determined in the calculation 
procedure. It was assumed that the mean density rise across the rotor 
was one half the stage rise. Then, at the rotor exit, for continuity 



Through the use of figure 2(b), the value of h2,h and therefore the 
value of h2 at any radiuB could be obtained for this rotor exit sta- 
tion. Actual and effective values of inlet -air angle and turning angle 
were calculated. 

The velocity vectors for the stator were determined by a similar 
procedure. For the Btator, however, there was a discrepancy introduced 
by the assumption that the velocity diagrams were symmetrical because 
the required change in diagram from stage to stage was not considered. 
Therefore, the stator diagram was so modified that the rotation removed 
at each radius in the stator Y was that required to provide the 
desired absolute rotational component X for the following rotor. The 
velocity diagram therefore deviated from symmetry to a limited extent. 
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Determination of pressure ratio and weight flow . - The total- 
pressure ratio of the compressor was established by consideration of the 
anticipated performance of inlet guide vanes and the exit vanes 
(designed to remove rotation remaining after the tenth stage) in addi- 
tion to the 10 stages for which static -pressure ratios were calculated. 



p 22 £22 £21 
P22 ?21 p l , 


£o £0 

P 0 Pi 


where 






r 

r-i 


P 22 

521 



+ 






n 

n-1 


For the inlet guide vanes, a polytropic expansion efficiency of 0.85 
was assumed; for the exit vanes, a polytropic compression efficiency of 
0.85 was used. Corresponding density relations were determined using 
these equations with the required changes in exponents. 
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Compressor -weight flow was established from 

W f = p-L h x U t A x 


where 



Since Wf , ho, and ao were unknown, the solution for pressure 
ratio and weight flow (assuming standard sea-level values for inlet 
stagnation conditions) required iteration which started with the assump- 
tion that 


and 


hQ = h x 


®0 * a l 


This provided a first approximation to p^_ and therefore Wf . Then 


p 0 


t,0 


1 - 0 +£ rfe) *o 2 


i 

ir-i 


^ Po^tAo 


and 



from which pi could be determined. A single iteration was found to 
give sufficiently accurate values of p-j_ and Wf . With the weight flow 

established, the conditions at station 22 were determined in similar 
manner. 
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AIRFOIL SELECTION" 

With, the values of stagger angle and turning angle established, 
compressor airfoils were selected for the various stages, incl uding 
the exit vanes. The design used the HACA 65-series airfoil scaled down 
to 10 percent thickness and modified by the addition of 0.0015 times 
chordal position to the thickness ordinate, designated in reference 7 
as the NACA 65-series blower blade section. 

In order to simplify construction, the blades were designed with 
constant camber and chord at all radii of a given blade row rather t han 
using a variable section blade. Thus the solidity varies inversely 
with the radius in any stage. 

In selecting blades, data from references 7, 8, and 9 were used. 
The guide vane data of reference 9 were used because no suitable data 
were available in the low inlet -air -angle range; interpolation was 
therefore necessary. The use of this zero inlet-air -angle data as an 
end point appeared to be reasonably accurate. 

The design charts of reference 8 provide a method of selecting 
optimum camber (represented by the free-stream lift coefficient for 
zero incidence Ci^i) and angle of attack a. The data were extended 
to cover low inlet-air angles and low solidities by a series of cross- 
plots, using the data of reference 9 as a guide. Reasonable consis- 
tency was obtained and the data were put in chart form (fig. 4) so 
that values of Cj^i and a could be selected for given values of 8, 
Ap, and a. 

With the desired solidity at the root and with a selected axial 
projection of each blade, an approximation to the number of blades in 
each row was made. The actual number of blades in each row was then 
determined from a consideration of blade exciting forces; a common 
den omina tor was avoided in adjacent rows. 

Through the use of figure 4, values of C L ± were determined at 

five radii in each blade row. It was found that the optimum blade sec- 
tions for this compressor had mini-nrim camber at the tip and maximum 
camber at the root. An intermediate value was chosen as the const ant - 
camber section of the blade row, thus obtaining optimum blading near 
the pitch radius and departing slightly at root and tip. The camber 
was chosen to maintain the same blade section for both rotor and 
stator wherever possible. 

With the blade section established, it was necessary to establish 
the angle of attack necessary to obtain the required turn ing at each 
radial station. Examination of the data of references 7, 8, and 9 
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indicated that these turn ing data could be roughly approx imat ed by a 
straight line variation of Aj3 with respect to a. The equation used 
was 


Ap = ma + b 

Charts (fig. 5) were prepared for the determination of m (as a 
function of p and c) and b (as a function of p, Cl, 1> and- cr) 
from the available data on 65-series blower blade sections. The maxi- 
mum difference between the cascade test data and the turning angle 
determined from the charts was of the order of 1°, with the nominal 
variation of the order of l/2°. With the use of these charts, the 
value of a required to accomplish the design turning at each of 
5 radial stations was determined. 


With the blade angle known and with the projected chord and number 
of blades specified, the actual chord and the solidity of each blade 
row could be determined. It was also possible to establish the actual 
value of C L from the relation 


c L 


2U t Y 

WmCf 


Reference 10 (table II ) gives the mean line data for the NACA 65- 
series blower blade sections with CL,i = 1.0. For the prescribed 
camber, the coordinates of the mean line and the tangent of the angle 
of the mean line were obtained by multiplying by the design value of 
CL,i* By the use of the data given in reference 10 (table I) for a 
10 percent thick HACA 65-series blower blade section, the coordinates 
of the airfoil were determined. 


Sheet-metal inlet guide vanes were selected to provide the turning 
as established by the inlet-stage velocity diagram. Blade sections were 
chosen by use of the methods subsequently published in references 11 
and 12, which considered both two-dimensional blade-element performance 
and secondary flow effects resulting from the radial distribution of 
design circulation. 


DESIGN CHARACTERISTICS 

The various design characteristics determined in the design pro- 
cedure are as follows (values given are for standard sea-level inlet 
conditions): 
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Over -all total-pressure ratio . 6.45 

Average stage total -pres sure ratio 1.205 

Weight flow, Ib/sec (no boundary layer allowance) 57.5 

Weight flow per unit frontal area, lb/sec 26.4 

Tip speed, ft/sec 869 

Tip speed, rpm 9959 

Over -all adiabatic efficiency 0.86 


The values of static -pres sure ratio per stage are plotted in fig- 
ure 6(a) . The discontinuities in the plot are a result of the change 
in the CjU limitation through the compressor. Ratio of static pres- 
sure ratio to the total pressure at the compressor inle t is shown in 
figure 6(b); total-pressure ratio is shown in figure 6(c). 

The average dimensionless axial velocity h and the average axial 
Mach number M entering each stage of the conpressor are shown in fig- 
tire 7. The distribution of axial velocity across the annulus is shown 
in figure 8. Figure 8(a) represents the dimensionless axial velocity 
at the entrance to each rotor; figure 8(b) represents the dimensionless 
axial velocity at the station corresponding to the stator entrance. 

The values of huh radius, blade chord, axial projection of the 
blade chord at the huh, number of blades, design camber, and maximum 
thickness of blades are given in table I for both rotor and stator 
rows. A summary of the velocity diagram data is given in table H for 
rotor and stator rows; values of Z, X, Y, h, W/U-fc, V/U-fc, W/a, V/a, 
and 0 are listed for five equally spaced stations at rotor and stator 
entrance. Changes through the given blade row (Ah and Ap) and blade 
data (a, Cpcr, and a) are also included. 


OVER-ALL TEST RESULTS 

The compressor was tested in accordance with the standard procedure 
of reference 13; characteristics as determined experimentally are pre- 
sented in figure 9. At design speed the peak compressor pressure ratio 
was 7.52, the maximum weight flow was 56.2 pounds per second, and the 
peak adiabatic efficiency was 0.83. The experimentally determined per- 
formance is discussed in greater detail in reference 14. 


SUMMARY OF DESIGN PROCEDURE 

The design of the 10 -stage axial -flow compressor was based on 
the use of a symmetrical velocity diagram with constant total enthalpy 
at all radii. The configuration chosen had a constant tip diameter of 
20 inches, an inlet hub-tip ratio of 0.55, and a hub diameter which 
increased at the rate required to maintain the selected limiting Mach 
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number for the blade rows. Aerodynamic limitations used included a 
relative Mach number limit of 0.7, a solidity of 1.1, and. a loading 
limit (CjjJ where is the lift coefficient and c is the blade- 

element solidity) which increased from 0.8 to 1.0 through the compressor. 
The axial velocity entering the compressor was chosen to approach 
closely maximum useful power input for the inlet stage . Axial -velocity 
distributions through the compressor were established by the simplified 
radial equilibrium approximation . Velocity diagrams were determined 
throughout the compressor, and 65-series blower blades were selected to 
accomplish the required turning from available cascade data. 

This design procedure with an assumed stage polytropic efficiency 
of 0.90 established the following design -point characteristics: a pres- 
sure ratio of 6.45 and an equivalent weight flow of 57.5 pounds per 
second at the design equivalent tip speed of 869 feet per second. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio 
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APPENDIX A 

DERIVATION OF POWER PARAMETER 

The useful power of the inlet stage of this design can "be approxi- 
mated "by the following equation, since the change in velocity across 
the stage is small: 


hp 




ni I 

Pi9i 

( p 3 

A p i 

J - 1.0 
' — 


550 


Transposing and dividing by 


Jtrt 2 a.i 


gives 


550 hp 


7^ p i 


itr. 


t“ a l 


Qi 

jtr t 2 ai 



1.0 


For given inlet conditions and compressor diameter. 


K(hp) = 


<*1 

2 

*r t a-L 



1.0 


where 


Ql 


(1 



3trt 2 ap 


Ut 

a l 



22 


KA.CA EM E52B18 


(for a given Clct limit, tlie average axial velocity h is a function 
of and Z^) and 


£3 

^1 


1.0 + (r-i) 



— 1 El 


n-1 


YZ 


as developed in step 11. 
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APPENDIX B 

DETERMINATION OF AXIAL VELOCITY 

The axial velocity at the various stations and radii through the 
compressor was determined from consideration of the simplified radial 
equilibrium relation 

dp p(xo t ) 2 

dr r r 


1 

K(p)r (xo t ) 2 

r 


Eor the symmetrical diagram (fig. 3(a)), at the rotor entrance 


XU+ 


_ ZU t YU t 


and since 


Therefore, 


Y = — 

z 


ZU + Y,U, 

XD t = _J i - _LI 

x 2 2Z 




Kp f 

Integrating hetveen the hub and a given radius yields 


/ r=i Izk 

“PH / 


^t 2 

Z 2 - Zh 2 

4 

2 


EYtloge^-lYt 2 ^-^ 
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But 


i 2 = ip r 


Therefore , 


12 r /a h \ 2 ~] Ut 2 
^1 1_ X " \~& ) ~ 4 


Ut 2 ( Z 2 - Zh 2 


ov t z 1 Yt " 

- 2Y t log e ^ - 5 p- 




?) 2 --(¥)(^) 2 ( 


Z 2 y loff J. It! . 

4 - Y t log e 2 ^ “ 4Z 2 + 4 , 


For the assumed condition of constant total enthalpy 


V h 2 + 2g - V 2 + 2g JCj. Rt 


X (h k 2 + X k 2) + _5_(Si) 2 . !j£ Ch2 + IS) + £ 


Substituting the sonic velocity relation gives 

/ Z 2 - Z >, 2 

t 2 + X 2 . h h 2 + X h 2 _ l _A_ _ Y t loge 


a _ hi + 1^!) 

Zh 4 z2 4 / 


With the substitution of Xfr = 


2 h - Y h 


, X = 




, and ZhXh = Yt; 


* - \\*t? + \ - T * Vh ^ 2J 


This equation established the distribution of axial velocity ahead of 
the rotor in terms of the axial velocity at the hub. 
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The value of Yfr for a symmetrical velocity diagram was. found "by 
the simplified lift consideration outlined in reference 4 as 


^ 

For design purposes, it was also necessary to establish the average 
dimensionless axial velocity for use in continuity relations. When no 
variation in density from hub to tip is assumed. 


h = 


1 - v 


''Zt 

Zh 


hZ dZ 


h = 


1 - Zj, 


r»i 


z 



+ 22^2 - 


2Z2 + 4Z h Y h log e 



dZ 


In order to expedite design, this equation was integrated graphi- 
cally and a design chart was prepared where h is shown as a function 
of Z^, h h , and Cja (fig. 2(a)). 

For the station at the stator entrance, the value of h was found 
similarly as 



%■ - Zh Y h lo ®e 


Z 

Zh 


A chart for the determination of h at the stator entrance was 
prepared for this station (fig. 2(b)). 
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TABLE I - BLADE DATA 


Stage 

Blade 

Hub 

radius 

at 

entrance 

(in.) 

Number 

of 

blades 

Blade 

chord 

(in.) 

Axial 
projection 
of chord 
at hub 
radius 

Design 

camber 

%i 

in « BSlSTOupl} 

I? 

Bl 

Rotor 

5.500 

25 

1.550 

1.506 

0.91 

10 

SI 

Stator 

5.822 

27 

1.550 

1.501 

.91 

10 

2 

Rotor 

6.131 

28 

1.510 

1.453 

0.93 

10 

2 

Stator 

6.420 

29 

1.510 

1.452 

.93 

10 

3 

Rotor 

6.696 

31 

1.465 

1.397 

0.95 

10 

3 

Stator 

6.952 

33 

1.465 

1.401 

1.025 

10 

4 

Rotor 

7.197 

Bl 

1.420 

1.348 

1.110 

10 

4 

Stator 

7.424 


1.420 

1.346 

1.110 

10 

5 

Rotor 

7.639 • 

— 

1.380 

1.300 

1.120 

10 

5 

Stator 

7.838 


1.380 

1.298 

1.120 

10 

6 

Rotor 

8.026 

41 

1.330 


1.120 

10 

6 

Stator 

8.199 

43 

1.330 


1.200 

10 

B 

Rotor 

8.361 

45 

1.280 

1.201 

1.280 

10 

H 

Stator 

8.508 

46 

1.280 

1.202 

1.280 

10 

8 

Rotor 

8.647 

47 

1.225 

1.149 

1.300 

10 

8 

Stator 

8.771 

49 


1.150 

1.300 

10 

9 

Rotor 

8.885 

51 

1.175 

1.101 

1.300 

10 

9 

Stator 

8.988 

53 

1.175 

1.103 

1.300 

10 

10 

Rotor 

9.085 

56 



1.300 

10 

10 

Stator 

9.174 

59 

fggg 


1.300 

10 


Exit vane 

9.262 

59 

1.120 

1.117 

0.660 

10 
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TABLE II - VELOCITY DIAGRAM DATA 



Station 

Entrance vector data 

Change 

through 

row 

- 

Blade data 


Z 

X 

Y 

h 

D t 

w 

a 

P' 

Ah 

A3 

a 

Cja 

a 

1 

0.5500 

0.1153 

0.3195 

0.7500 

0.8669 

0.7000 

30.10 

0.0980 

22.36 

16.44 

0.7562 

1.122 


.6625 

.1987 

.2652 

.7259 

.8615 

.6952 

32.58 

.0604 

18.40 

13.46 

.6425 

.931 


.7750 

.2742 

.2267 

.6883 

.8513 

.6862 

36.04 

.0257 

15.03 

11.24 

.5659 

.796 


.8875 

.3448 

.1980 

.6357 

.8359 

.6726 

40.49 

-.0092 

11.66 

9.20 

.5132 

.697 


1.0000 

.4122 

.1757 

.5646 

.8151 

.6545 

46.20 

-.0494 

7.54 

6.80 

.4774 

.617 

3 

0.6131 

0.1449 

0.3233 

0.7479 

0.8824 

0.6969 

32.05 

0.0925 

22.27 

16.34 

0.7595 

1.098 


.7098 

.2153 

.2792 

.7242 

.8769 

.6922 

34.33 

.0589 

18.96 

13.86 

.6703 

.948 


.8065 

.2804 

.2458 

.6902 

.8679 

.6843 

37.32 

.0272 

15.97 

11.76 

.6061 

.835 


.9032 

.3419 

.2194 

.6451 

.8552 

.6735 

41.03 

-.0052 

12.91 

9.91 

.5588 

.746 


1.0000 

.4009 

.1982 

.5868 

.8386 

.6592 

45.59 

-.0415 

9.27 

7.79 

.5248 

.673 

5 

0.6696 

0.1700 

0.3297 

0.7533 

0.9040 

0.6971 

33.56 

0.0829 

22.07 

16.05 

0.7645 

1.079 


.7522 

.2294 

.2935 

.7311 

.8988 

.6926 

35.57 

.0530 

19.26 

14.02 

.6939 

.961 


.8348 

.2852 

.2645 

.7014 

.8912 

.6863 

38.09 

.0241 

16.63 

12.28 

.6400 

.866 


. .9174 

.3384 

.2407 

.6636 

.8808 

.6775 

41.11 

-.0055 

13.90 

10.60 

.5983 

.788 


1.0000 

.3896 

.2208 

.6167 

.8677 

.6666 

44.71 

-.0381 

10.76 

8.85 

.5667 

.723 

7 

0.7197 

0.1709 

0.3779 

0.7591 

0.9367 

0.7042 

35.87 

0.0770 

24.32 

17.51 

0.8637 

1.067 


.7898 

.2227 

.3444 

.7407 

.9328 

.7010 

37.44 

.0473 

21.66 

15.59 

.8006 

.972 


.8599 

.2718 

.3163 

.7169 

.9273 

.6965 

39.36 

.0179 

19.06 

14.00 

.7498 

.894 


.9300 

.3188 

.2925 

.6873 

.9198 

.6903 

41.65 

-.0122 

16.57 

12.31 

.7093 

.826 


1.0000 

.3640 

.2720 

.6517 

.9106 

.6827 

44.30 

-.0446 

13.35 

10.59 

.6768 

.768 

9 

0.7639 

0.1899 

0.3842 

0.7636 

0.9554 

0.6973 

36.94 

0.0715 

24.13 

17.33 

0.8673 

1.064 


.8230 

.2332 

.3566 

.7471 

.9518 

.6943 

38.29 

.0457 

21.90 

15.69 

.8169 

.988 


.8820 

.2746 

.3328 

.7267 

.9471 

.6906 

39.89 

.0203 

19.71 

14.28 

.7750 

.922 


.9410 

.3146 

.3119 

.7024 

.9412 

.6859 

41.73 

-.0058 

17.43 

12.91 

.7399 

.864 


1.0000 

.3533 

.2935 

.6737 

.9340 

.6801 

43.83 

-.0333 

14.94 

11.49 

.7109 

.813 

11 

0.8026 

0.2045 

0.3937 

0.7775 

0.9810 

0.6946 

37.57 

0.0645 

23.92 

17.08 

0.8711 

1.081 


.8520 

.2406 

.3709 

.7632 

.9780 

.6923 

38.70 

.0425 

22.07 

15.73 

.8309 

1.019 


.9013 

.2754 

.3506 

.7464 

.9742 

.6893 

39.99 

.0205 

20.24 

14.58 

.7961 

.963 


.9507 

.3092 

.3324 

.7268 

.96 95 

.6857 

41.44 

-.0018 

18.34 

13.44 

.7662 

.914 


1.0000 

.3420 

.3160 

.7043 

.9639 

.6815 

43.05 

-.0250 

16.33 

12.32 

.7405 

.868 

13 

0.8361 

0.1935 

0.4492 

0.7952 

1.0224 

0.7021 

38.97 

0.0592 

26.21 

18.69 

0.9716 

1.094 


.8771 

.2245 

.4282 

.7843 

1.0203 

.7005 

39.77 

.0384 

24.51 

17.45 

.9354 

1.045 


.9181 

.2545 

.4091 

.7716 

1.0177 

.6985 

40.70 

.0176 

22.83 

16.38 

.9037 

.999 


.9591 

.2838 

.3916 

.7572 

1.0147 

.6963 

41.73 

-.0035 

21.10 

15.30 

.8752 

.956 


1.0000 

.3122 

.3756 

.7409 

1.0110 

.6935 

42.87 

-.0252 

19.30 

14.29 

.8504 

.917 

15 

0.8647 

0.2022 

0.4604 

0.8129 

1.0487 

0.6960 

39.18 

0.0550 

26.07 

18.87 

0.9743 

1.060 


.8986 

.2278 

.4430 

.8039 

1.0470 

.6948 

39.84 

.0376 

24.69 

17.88 

.9453 

1.020 


.9324 

.2527 

.4270 

.7937 

1.0450 

.6934 

40.58 

.0202 

23.33 

17.06 

.9192 

.983 


.9662 

.2771 

.4120 

.7824 

1.0426 

.6916 

41.37 

.0025 

21.92 

16.21 

.8951 

.949 


1.0000 

.3010 

.3981 

.7698 

1.0399 

.6897 

42.24 

-.0153 

20.49 

15.30 

.8734 

.916 

17 

0.8885 

0.2074 

0.4737 

0.8366 

1.0788 

0.6922 

39.15 

0.0504 

25.99 

18.72 

0.9772 

1.074 


.9164 

.2286 

.4593 

.8293 

1.0775 

.6913 

39.68 

.0358 

24.88 

17.92 

.9537 

1.041 


.9443 

.2493 

.4457 

.8212 

1.0758 

.6901 

40.24 

.0213 

23.76 

17.18 

.9319 

1.010 


.9722 

.2697 

.4329 

.8123 

1.0740 

.6888 

40.86 

.0067 

22.63 

16.45 

.9118 

.982 


1.0000 

.2896 

.4209 

.8027 

1.0720 

.6874 

41.51 

-.0080 

21.49 

15.85 

.8933 

.954 

19 

0.9085 

0.2087 

0.4912 

0.8710 

1.1174 

0.6933 

38.78 

0.0498 

26.01 

18.49 

0.9780 

1.099 


.9314 

.2261 

.4792 

.8652 

1.1162 

.6925 

39.19 

.0379 

25.13 

17.95 

.9590 

1.072 


.9543 

.2433 

.4677 

.8590 

1.1151 

.6918 

39.61 

.0260 

24.24 

17.30 

.9410 

1.047 


.9772 

.2603 

.4567 

.8522 

1.1137 

.6908 

40.08 

.0141 

23.36 

16.67 

.9239 

1.022 


1.0000 

.2769 

.4463 

.8450 

1.1122 

.6898 

40.56 

.0022 

22.46 

16.10 

.9081 

.998 
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TABLE II - VELOCITY DIAGRAM DATA - Concluded 



Station 

Entrance vector data 

Change 

through 

row 

Blade data 


Z 

X 

Y 

h 

7 

V 

a 

8 

Ah. 

Af3 

a 

CjO 

a 

2 

0.5822 

0.1209 

0.3211 

0.8312 

0.9414 

0.7578 

28.00 

-0 .0779 

18.88 

13.74 

0.7639 

1.143 


.6866 

.1990 

.2722 

.7718 

.9043 

.7247 

31.40 

-.0410 

16.17 

11.91 

.6617 

.970 


.7911 

.2703 

.2363 

.7026 

.8661 

.6912 

35.79 

-.0063 

14.57 

10.87 

.5906 

.842 


.8955 

.3371 

.2087 

.6194 

.8256 

.6559 

41.39 

.0298 

13.95 

10.37 

.5401 

.744 


1.0000 

.4009 

.1869 

.5152 

.7816 

.6183 

48.77 

.0716 

14.43 

10.41 

.5049 

.666 

4 

0.6420 

0.1491 

0.3262 

0.6240 

0.9513 

0.7476 

29.98 

-0.0649 

18.87 

14.04 

0.7667 

1.086 


.7315 

.2149 

.2863 

.7692 

.9181 

.7189 

33.09 

-.0318 

16.84 

12.67 

.6865 

.952 


.8210 

.2761 

.2551 

.7067 

.8841 
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(a) At rotor entrance. 

Figure 2. - Chart for determination of average dimensionless axial 

velocity. 
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Hub-tip ratio, Z h 
(b) At stator entrance. 

Figure 2, - Concluded. Chart for determination of average dimensionless 

axial veloolty. 
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Blade -element solidity, a 


Figure 5. 


( a ) Determination of m in equation Aft — wy, + t • 

- Design chart for HACA 65-series hlcwer "blade sections 
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Figure 6. - Concluded. Prea sure -ratio variation through 
compressor. 
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Figure 7. - Variation of average dimensionless axial velocity 
and average axial Mach, number through compressor. 
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(a) At rotor entrance. 
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Figure 8. - Distribution of dimensionless axial velocity. 
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Figure 9. - Over-all performance characteristics of 10-stage axial -flow 

research compressor. 
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